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Abstract 

We present some exact results on percolation properties of the Ising model, 
when the range of the percolating bonds is larger than nearest-neighbors. We 
show that for a percolation range to next-nearest neighbors the percolation 
threshold T p is still equal to the Ising critical temperature T c , and present 
the phase diagram for this type of percolation. In addition, we present Monte 
Carlo calculations of the finite size behavior of the correlated resistor network 
defined on the Ising model. The thermal exponent t of the conductivity that 
follows from it is found to be t = 0.2000 ± 0.0007. We observe no corrections 
to scaling in its finite size behavior. 

PACS numbers: 64.60.Ak, 64.60.Fr, 05.50.+q 
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I. INTRODUCTION 



The connection between percolation and the Ising model has been a popular subject for 
a long time. One considers so-called Ising clusters made up of nearest-neighbor spins with 
the same spin value. The connectivity behavior of these clusters is called correlated site 
percolation, as the probability distribution of the percolating and non-percolating sites is a 
correlated one. 

The interest in this problem arose because these Ising clusters were believed to have 
the same properties as the droplets in the droplet model |]], i.e. they should diverge at the 
Ising critical point with the same critical exponents as those of the Ising model. It became 
clear that they did indeed diverge [[| at the Ising critical point, but not || with Ising 
exponents. An alternative cluster definition was needed to have clusters with the properties 
of droplets in the droplet model. These clusters, called Coniglio-Klein clusters are defined 
by putting bonds between each pair of nearest neighbor up-spins, but now with a probability 
p = 1 — exp(— 2K), where K is the Ising coupling. Not all bonds of the Ising clusters appear 
in the Coniglio-Klein clusters, such that the latter are, in that sense, 'smaller' than Ising 
clusters. The Coniglio-Klein clusters display the right critical behavior: their diverge 
at the Ising critical point, their linear size diverges as the Ising correlation length, and the 
mean cluster size behaves as the susceptibility. 

Both the Ising clusters and the Coniglio-Klein clusters have their percolation point at 
the Ising critical temperature, albeit with different critical behavior. The full picture of this 
cluster behavior emerged |5]|| when the behavior of both types of clusters was identified 
with the phase diagram of the g-state dilute Potts model in the limit q — > 1. The tricritical 
point in this phase diagram describes the behavior of the Ising clusters. This tricritical point 
falls in the same universality class as the Ising critical point [|7| in the sense that the central 
charge is c = 1/2, but the critical exponents involved in the behavior of the Ising clusters 
do not fit into the Kac-table ||; they correspond to half- integer values of the unitary grid. 
The Coniglio-Klein clusters are described by the 1 + 1 state symmetric fixed point in this 
phase diagram. 

Our motivation to reconsider the problem of correlated site percolation is not the droplet 
theory of the Ising model but arises from a study of correlated resistor networks; see further 
below. In this paper, we consider Ising clusters that are made up of bonds with a larger 
percolation range, that is, bonds are placed between nearest neighbor up spins, but also 
between next- and further-neighboring pairs of equal spins. Let us, throughout this paper, 
denote the clusters consisting of bonds between nearest and next-nearest neighbor spin pairs 
by nnn-clusters. In this language the Ising clusters are nn-clusters. 

It is immediately clear that, if one considers bonds with a longer and longer percolation 
range, clusters get bigger and the percolation threshold eventually will move to a temperature 
T p that is lower than T c , the Ising critical temperature. In that case, the type of correlation is 
expected to be random percolation. In the limit of percolating bonds with an infinite range, 
the percolation temperature moves to T p = and there is crossover to classical critical 
behavior ||. In three dimensions this effect of a shift in the percolation threshold occurs 
already with Ising clusters (nn-clusters): the percolation threshold lies at a temperature a few 
percent below T c |]IU|| , whereas the Coniglio-Klein clusters have their percolation threshold 
at T c . In two dimensions it is known, as stated above, that for nn-clusters T p coincides with 
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T c but it was believed JTTJ that already for nnn-clusters T p < T c . However, in this paper we 
shall show that T p = T c for nnn-clusters. 

These alternative cluster definitions can be useful in some applications of correlated per- 
colation. In another paper fl2|| , we present a model, based on correlated percolation, to 
explain the experimental results for colossal magnetoresistance. The latter phenomenon is 
presently a hot topic in solid state physics ||13|| . Our model is a correlated resistor net- 
work, obtained by replacing bonds with resistances yielding an effective resistance as an 
Ising expectation value. In the present work, we present the technical analysis of the cor- 
related percolation model with percolating bonds having a longer range of interaction. In 
particular, the resulting phase diagram is vital for understanding the experimental results 
of colossal magnetoresistance. The correlated resistor network has, to our knowledge, never 
been studied in the literature. We performed Monte Carlo calculations to measure the 
critical exponents of the CRN. Also these calculations are presented in this paper. 



II. THE MODEL 

We will be concerned with the usual Ising model on a square lattice with Hamiltonian 

H = -K^2s i s j -hY;s j , 

(ij) 3 

where K is the inverse temperature and h is the magnetic field. The first summation is 
over nearest neighbors only. Note that we are going to define percolating bonds that have a 
longer range than nearest neighbors only; the Ising Hamiltonian however, interacts only via 
nearest neighbors. 

First consider the nn-clusters by putting bonds between all neighboring pairs of up spins, 
such as the clusters made of the black spins in figure [I]. In the same figure, also the nnn- 
clusters are illustrated but now for the white spins. Here bonds are put between next-nearest 
neighboring pairs of white spins as well. 

The percolation phase diagram for the nn-clusters is known || and shown in figure |2|. 
The thick solid line is a critical percolation transition; the end point at T = oo corresponds 
to random percolation and lies at the value of h corresponding to the percolation threshold 
|TJ|] p c pa 0.5927 for random percolation. This value denotes the density of up spins. The 



corresponding value of h is h = 0.188. The critical percolation line is in the universality class 
of random percolation, described by the critical q = 1 state Potts model. The line merges 
smoothly with the T-axis at the Ising critical point. For percolation, this point turns out 
to be a tricritical point; it is the tricritical point of the q = 1 state dilute Potts model, 
where, apart from the usual Potts spins, also vacancies are allowed. 

From this phase diagram the corresponding diagram for nnn-clusters can be derived. 
Figure |l] illustrates the theorem [ljj that will be needed for this derivation. It can easily be 



seen from this figure that there exists a geometrical relation between the black nn-clusters 
and the white nnn-clusters. The theorem states that every face of a black nn- cluster is 
either empty or is wholly occupied by one and only one white nnn-cluster. A face is a 
closed area surrounded by elementary loops consisting of the bonds of one cluster. Let us 
introduce the following notation: let fs, cb, &b and % be the number of faces, clusters, 
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bonds and sites respectively that correspond to the black spins and nn-clusters in a certain 
spin configuration. For the white spins and nnn-clusters, these quantities are f^, c^, b^ 
and n^y respectively, where the star denotes the fact that it concerns nnn-clusters. The 
theorem then states that 

/b = f B ^ + C Wi 

where denotes the number of empty black faces. The number of faces, bonds, sites and 
clusters is furthermore coupled via Euler's equation 

Ib = b B - n B + c B + 1. 

This relation is easily derived by induction and holds for non-cyclic boundary conditions, 
e.g., all spins on the boundary are black. Combining these two relations yields 

c* w = c B + b B - n B - f B 0) + 1. 

Thus this relation expresses the number of white nnn-clusters in terms of the number 
of black nn-clusters c B . Apart from the numbers of clusters it involves only locally defined 
quantities: bonds, sites and empty faces. 

The above relations are purely geometrical and hence are completely independent of 
the probability distribution of the black and white spins. Therefore they are valid as Ising 
expectation values as well: 

(c* w ) = (cb) + (b B ) - (n B ) - (/£>) + 1, (1) 

where the brackets denote the expectation values. Note that the numbers of bonds, sites 
and empty black faces are simply local Ising operators, that is 

<&£> = ^E<( 1 + ^)( 1 + ^)>> 

{ij) 

which amounts to counting the number of black bonds. In the same way 

3 

counts the number of black spins and 

Ub) = yq E ^ + + s i)0- + Sk)0- + Si)), 
MM) 

counts the number of empty, black faces, where (ijkl) denotes a summation over each 
elementary plaquette of the square lattice. 

The expectation value of the number of clusters plays the role of the free energy in 
a percolation problem ||15|| . It becomes non-analytic at a percolation transition. From 
equation (|I|) we see that C3jI1 only become critical when c B is critical, that is, when the 
black spins are at their percolation threshold, or when the Ising expectation values become 
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non-analytic, that is, at the Ising critical point and at the coexistence line T < T c , h = 0. 
This immediately yields the phase diagram for percolation of the black nnn-clusters: it is 
the mirror image of that of the nn-clusters with h replaced by —h. This phase diagram is 
shown in figure |[ 

It is somewhat surprising that extending the range of percolating bonds to next-nearest 
neighbors does not alter the percolation threshold at the T-axis. It was expected [llj that 
for larger ranges than nearest neighbor percolating bonds the percolation temperature T p 
would become less than T c . Due to the above geometrical relations this is not the case. 
For even larger percolation ranges, however, there are no such relations and we expect the 
percolation temperature T p to drop below the Ising critical temperature T c . 

The topology of the phase diagram then changes. There is no tricritical point anymore, 
and the critical percolation line will end somewhere at the T-axis below T c . Beyond this 
point, there still is a first order transition for percolation. We expect the percolation point 
T p in that case to be a critical endpoint. Such a point is expected |R| to have in addition to 
the critical exponents of the universality class of the critical line also a first order exponent 
y = 2. Indeed, when the range of percolating bonds becomes very large, it will eventually 
become larger than the Ising correlation length. In such a case, correlations in the probability 
distribution of empty and occupied sites are only present on a much smaller scale than the 
range of percolating bonds. This strongly suggests that the universality class of the endpoint 
of the critical line is that of random percolation. 

Eventually, when the percolation range approaches infinity, percolation is believed to 
exhibit classical critical behavior HI. 



III. SCALING ANALYSIS 

From universality, we expect both types of critical behavior, nn-percolation and nnn- 
percolation, to be in the same universality class. In the light of equation ([]]), this statement 
is less obvious than it seems. The singular behavior of the 'free energy' c? w of nnn-clusters 
is expressed in cb but also in Ising operators. Hence in addition to the critical behavior of 
cb also Ising exponents show up. In particular, when is considered as a function of the 
scaling field Ui, in addition to the exponent 2/y that is expected to describe the non-analytic 
behavior of the free energy eg, also the magnetization exponent 1/8 is present. 

Nn-percolation is in the universality class of the dilute g-state Potts model in the limit 
of q — > 1. Its Hamiltonian HQ is 

H = -L^ n i n i ~ A ^2 n i ~ J n i n i( 5 ^^ ~ 1 ) ~ ^'('W ~ X )- ( 2 ) 

(ij) j (ij) 3 

The variables rii = 0,1 are the Potts lattice gas variables. Potts spins Oi = 0, . . . , q are 
present on the sites where = 1. 

For q — 1 the Hamiltonian becomes equal to the Ising Hamiltonian and completely 
independent of J and H. Substituting Si = 2rii — 1 turns the lattice gas variables into Ising 
variables. The Hamiltonian then becomes, apart from a constant, 

H = -K y £s i s j -hY;s j , 

(ij) j 
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with K = L/4 and h = 2L+A/2. This means that the free energy resulting from equation (0) 
for q = 1 is equal to that of the Ising model. The full free energy is 

f(L,A,J,H,q) = - lim -|- In (L, A, J, H, q) 

N^>oo iv 



where iV denotes the number of sites on the lattice. The generating function [17] for perco- 
lation is 



c(L,A,J,H) = ^ L ' A - J - H ^ 



dq 



9=1 



This quantity c is, with H = 0, the expectation value of the number of clusters. It plays the 
role of the free energy in percolation problems. The limit J — > oo corresponds to nn-clusters, 
and J = 2K yields the Coniglio-Klein clusters. 

In the language of the renormalization group, the tricritical point in the q = 1 dilute 
Potts model has four relevant exponents and corresponding scaling fields, two thermal and 
two magnetic ones. Expressing the free energy in terms of these scalings fields Ui, the scaling 
relation is 

f( Ul , ...,u 4 ,q)= b- 2 f(b y i Ul , . . . , 6 W « 4 , q), (3) 

where b is the renormalization length, and the tricritical point is located at u\ = ■ ■ ■ = 1x4 = 0. 
Note that the 'field' q, which is a symmetry parameter, cannot change under renormalization. 
Differentiating this free energy with respect to q yields the percolation free energy. Applying 
this to the above scaling relation yields, apart from the direct derivative to q, also derivatives 
with respect to the thermal scaling fields u\ and U2 which are 

u x = K- K c (q) + ■■■ (4) 
u 2 = h- h c (q) H (5) 

when J = 00. The location of the tricritical point thus depends on the value of q. The 
remaining two fields M3 and M4 are magnetic scaling fields; they correspond to Potts-like 
magnetic fields. The Potts-model can only be critical if those fields are zero, regardless of 
the value of q, which means that the derivative of these fields with respect to q yields zero. 

From these remarks it follows that in the expression for c several derivatives are present: 
a direct derivative with respect to q and derivatives with respect to the thermal scaling 
fields: 



9f _|_ df dui df du 2 



dq du\ dq du 2 dq 



(6) 



q=l 



Applying equation (0) to this expression and taking the limit of q — > 1 yields the critical 
behavior of c. Note that the two last terms of the right hand side of (||) are with q — 1 
just derivatives of the Ising free energy: they give the Ising energy and magnetization, and 
thus yield Ising critical behavior. The values of the exponents y\ and y 2 are the Ising values 
Hi = 1 and y 2 = 15/8. This gives for the singular behavior of c as a function of the thermal 
fields U\ 
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(7) 



The most relevant of these exponents is (2 — 2/2 ) = 1/8) the exponent of the Ising order 
parameter. The singular behavior of c thus is described by an exponent 1/8, contrary to 
what is expected for the behavior of a free energy-like quantity: the value of y\ — 1 suggests 
a critical exponent 1jy\ = 1 for the percolation free energy c. This exponent indeed occurs 
in the expression for c but is dominated by the exponent 1/8. This is due to the fact that 
c is not a true free energy but a derivative of a free energy with respect to a symmetry 
parameter. 

The order parameter P (the density of sites in the infinite cluster) and the 'susceptibility' 
(the mean cluster size) arise from differentiating c with respect to H. The free energy / 
itself becomes completely independent of H in the limit for q — > 1, which means that only 
the derivatives of df/dq with respect to H yield a non-zero result. The other terms in 
equation (||) vanish upon taking the derivative. This means that the critical behavior of the 
order parameter and the susceptibility is not affected by the 'wrong' critical behavior of c. 

From the scaling results (|7|) of c we can derive the critical behavior of the nnn-clusters 
using equation ([I]). The dominating exponent in the right hand side of (^ is again 1/8, which 
appears in c but also in the Ising magnetization (ris)- This shows that the critical behavior 
of c and c*, that is, of nn-clusters and nnn-clusters, is governed by the same exponents. 

The same must be true for the order parameter P and the susceptibility of nnn- 
percolation. Again these quantities arise from differentiating the free energy c* with respect 
to the Potts field H . Although c* is different from c in its dependence on H, both amount 
to a dilute q = 1 states Potts model. With the identification of the critical behavior of c 
and c* it follows that at T c both free energies must be at a tricritical point in the full phase 
diagram of this model. 

We conclude that percolation of nn-clusters and of nnn-clusters is in the same universality 
class. The critical behavior of the percolation quantities is governed by the same set of critical 
exponents. 



IV. THE CORRELATED RESISTOR NETWORK 



Relevant for our work on colossal magnetoresistance is the behavior of the Ising 
model as a resistor network. A percolation model is turned into a resistor network [18| 
by replacing the bonds with resistors. Non-percolating bonds are left empty (that is, have 
infinite resistance), bonds that are present get a unit resistance. This can be done both for 
nn-clusters and nnn-clusters. The assignment of resistors is depicted in figure [|. 

For random percolation the resistor network is a random resistor network. The corre- 
sponding correlated resistor network has as far as we know never been studied. In this 
section we present our calculations on the correlated resistor network. The interest in re- 
sistor network problems is in the expectation value of the overall resistance of the (infinite) 
lattice. To be more precise: consider a lattice consisting of L x L spins, where L eventually 
is sent to infinity, and keep the lower row at a fixed potential V = and the upper row at 
V = 1. The interest is in the overall conductance a of the lattice, which is in this case equal 
to the expectation value of the current. 
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The phase diagram of the resistor network is of course the same as that of its percolation 
counterpart: when there is percolation, the conductance is finite, and the conductance is 
zero when there is no percolation. Experimental results for colossal magnetoresistance show 
a sharp increase in the resistance as a function of temperature at or nearby the Curie point 
T c . The resistance drops sharply both above and below this point. From our phase diagrams 
it is directly clear that the diagram of figure |2| is ruled out. So conduction via next-nearest 
neighbor bonds should at least be present to produce the right phase diagram. 

Turning, however, to critical exponents the situation is different. In this case it is the 
exponent t governing the vanishing conductance a upon approaching the percolation thresh- 
old: 

o(T) ~ \T - T p \ l for T -> T p . 

Based on universality (again confirmed for the percolation exponents) one may well assume 
that the exponent t is the same for nn- and nnn-networks. Hence we studied the (simpler) 
case of the resistor network with only nn-clusters. 

The random resistor network is a notoriously unsolved problem in statistical physics, 
but there are good numerical estimates of the exponent t. The best estimate [[19| in two 



dimensions known to us is t = 1.299 ± 0.002. To obtain the value of t for the correlated 
resistor network in the Ising case, at the tricritical percolation point, we performed Monte 
Carlo calculations at the Ising critical point, and calculated the Ising expectation value of 
the conductance for different system sizes. We used the Wolff-algorithm [2IJ for the Monte 



Carlo part, and the multigrid method of Edwards, Goodman and Sokal [£1[], based on the 



standard code amg1r4 | 2"^fl , to calculate the conductance of a spin configuration. 

To test our program, we calculated the exactly known @J^] exponent of the order pa- 
rameter. The order parameter of percolation is the density of sites P in the infinite cluster. 
From finite size scaling, it follows that this quantity scales with the linear system size L as 

P(L) ~ L- 2+Vh for L -> oo, (8) 

where yh is the most relevant magnetic eigenvalue. In our Monte Carlo runs, we measured 
the number of sites in the 'spanning cluster', the cluster that extends over the lattice and 
thus allows for conductance. With the scaling equation, the behavior of P as a function of 
the system size yields an estimate of y^. We calculated P with the system size L running 
from 7 up to 350, and the data are plotted in figure |5|. The figure, a log-log plot, shows that 
the system sizes are too small to exhibit the behavior of equation (|8]); corrections to scaling 
have to be included. We did this, and found yh = 1.945 ± 0.005 and a correction to scaling 
exponent uj = 0.96 ± 0.08. The exact result is yh = 187/96 ~ 1.948. Our estimate thus 
agrees well within the error bars. 

It is believed that a similar finite size scaling relation is valid for the conductance a. It 
should scale with the linear system size L as 

a{L) ~ IT*/", (9) 

where v is the percolation exponent of the correlation length. We calculated the conductance 
again with L running from 7 up to 350. The data are plotted in figure || the log-log plot 
almost shows a straight line. In fitting the data to equation ([|), we tried to include a 
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correction to scaling term, but, due to the almost perfect scaling behavior, this did not yield 
sensible results. Therefore we performed the fit against equation (§) without additional 
terms, yielding the value t/u = 0.2000 ± 0.0007. Due to the absence of the correction to 
scaling term the error in this result might be an underestimation of the actual error. 

The exponent v of percolation at the Ising critical point is v — 1 for the direction parallel 
to the T-axis and v = 8/15 for the other directions. That means that the exponent t that 
governs the vanishing conductance at T c is t ~ 0.200 for the temperature direction and 
t m 0.107 for the field direction. This is a surprisingly low result, as compared to the t value 
of the random resistor network, t ~ 1.30. The presence of critical correlations thus strongly 
influences the value of the conductance exponents. 

Equation @ relies on the validity of finite size scaling, and it is in this case not a 
priori clear that it is valid, since there does not exist a rigorously defined renormalization 
transformation for the conductance. Therefore we checked the consistency of the value of 
the exponent by calculating the behavior of the resistance in the neighborhood of T c . The 
conductance then scales as o ~ \T — T c |*. The value of t obtained from this behavior was 
indeed seen to extrapolate to the value 0.20. We believe therefore that equation (|9|) must 
be valid. 

ACKNOWLEDGMENTS 

We thank Alan Sokal for providing us with the code for the resistance calculations, which 
saved us a lot of work; Bernard Nienhuis for discussions on correlated percolation; and Erik 
Luijten for discussions on the Monte Carlo part. 



9 



REFERENCES 



* Electronic address: paulbOtvs . kun . nl 
[1] Fisher M E 1967 Physics 3 255 

[2] Coniglio A, Rosanna Nappi C, Peruggi F and Russo L 1977 J. Phys. A 10 205 
[3] Sykes M F and Gaunt D S 1976 J. Phys. A 9 2131 
[4] Coniglio A and Klein W 1980 J. Phys. A 13 2775 
[5] Murata K K 1979 J. Phys. A 12 81 

[6] Stella A L and Vanderzande C 1989 Phys. Rev. Lett. 62 1067 

[7] Nienhuis B 1987 in Phase Transitions and Critical Phenomena vol 11, ed. C Domb and 

J L Lebowitz (New York: Academic) p 1 
[8] Cardy J L 1987 in Phase Transitions and Critical Phenomena vol 11, ed. C Domb and 

J L Lebowitz (New York: Academic) p 55 
[9] Gouker M and Family F 1983 Phys. Rev. B 28 1449 
[10] Miiller-Krumbhaar H 1974 Phys. Lett. 50A 27 
[11] Jan N, Coniglio A and Stauffer D 1982 J. Phys. A 15 L699 
[12] Bastiaansen P and Knops H (unpublished). 

[13] See, e.g., Kusters R M, Singleton J, Keen D A, McGreevy R and Hayes W 1989 Physica 
B 155 362; Chahara K, Ohno T, Kasai M and Kozono Y 1993 Appl. Phys. Lett. 63 
1990; McCormack M, Jin S, Tiefel T H, Fleming R M, Phillips J M and Ramesh R 
1994 Appl. Phys. Lett. 64 3045; Schiffer P, Ramirez A P, Bao W and Cheong S-W 1995 
Phys. Rev. Lett. 75 3336; Shimakawa Y, Kubo Y and Manako T 1996 Nature 379 53; 
and references herein. 

[14] Ziff R M and Sapoval B 1986 J. Phys. A 19 LI 169 

[15] Essam J W 1972 in Phase Transitions and Critical Phenomena vol 2, ed. C Domb and 

M S Green (New York: Academic) p 218 
[16] Berker A N and Wortis M 1976 Phys. Rev. B 14 4946 
[17] Fortuin C M and Kasteleyn P W 1972 Physica 57 536 
[18] See, e.g., Kirkpatrick S 1973 Rev. Mod. Phys. 45 574 
[19] Normand J-M, Herrmann H J and Hajjar M 1988 J. Stat. Phys. 52 441 
[20] Wolff U 1989 Phys. Rev. Lett. 62 361 

[21] Edwards R G, Goodman J and Sokal A D 1988 Phys. Rev. Lett. 61 1333; Goodman J 

and Sokal A D 1989 Phys. Rev. D 40 2035 
[22] Ruge J and Stiiben K 1987 in Multigrid Methods, ed. S F McCormick (Philadelphia: 

SIAM) 



10 



FIGURES 




FIG. 1. Illustration of the theorem presented in the text. On the sites of a lattice, black and 
white spins are placed. Spins on the boundary are always black. Bonds between black spins are 
present between each pair of black spins that are nearest neighbors. The same goes for the white 
spins, but here bonds are drawn if the spins are next-nearest neighbors as well. It can easily be 
inferred that a face of a black cluster is either empty or wholly occupied with one and only one 
white cluster. Notice that the inclusion of next-nearest neighbor bonds is essential. 
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non-percolating 



FIG. 2. The phase diagram for percolation of Ising clusters. Up spins are considered percolat- 
ing, and clusters are defined by putting bonds between each neighboring pair of up-spins. T c is the 
Ising critical point, and T and h are temperature and magnetic field, respectively. The thick, solid 
line is a critical line of percolation that is in the universality class of random percolation. T c is a 
tricritical point for percolation, and the dotted line is a first-order transition, for Ising as well as 
for percolation. 
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non-percolating 



FIG. 3. The same phase diagram as in figure ||, but now for nnn-clusters. Percolating bonds 
are drawn between pairs of up spins that are nearest or next-nearest neighbors. The phase diagram 
is the same as that for nn-clusters but with h replaced with —h. 
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(a) (b) 

FIG. 4. The assignment of resistors to a spin configuration. Assignment for nn-clusters is 
shown in (a), for nnn-clusters in (b). Here black spins (spin up) are considered percolating, so 
bonds between two black spins have a unit resistance. Bonds between black and white spins or 
between two white spins get an infinite resistance. 
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FIG. 5. Log-log plot of the density of sites P in the spanning cluster at the Ising critical point, 
as a function of the linear system size L. The plot shows a far from straight line, meaning that 
corrections to scaling are important. The dashed line is the result of our fit against the function 
P(L) = aiL- Q (l + a 2 L- /3 ). The values of a and are a = 0.055 ± 0.005 and (3 = 0.92 ± 0.07. 
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FIG. 6. Log-log plot of the conductance a of the lattice at the Ising critical point, as a function 
of the linear system size L. The plot shows almost a straight line. Fitting the scaling behavior 
with a correction to scaling term did not yield sensible results. The dashed line is the result of our 
fit against the function a(L) = aL~ a , giving a value of a = 0.2000 ± 0.0007. 
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